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Abstract.  This  report  studies  the  coordinated  tracking  problem  where  a  group  of  followers 
intercepts  a  dynamic  leader.  It  is  shown  in  this  paper  that  reduction  of  inter-agent  communication 
is  obtained  and  improved  performance  is  achieved  by  implementing  dynamical  models  of  the  leader 
and  followers  and  by  using  an  event-triggered  control  strategy  that  requires  each  agent  to  send 
measurement  updates  only  when  necessary.  Communication  delays  are  considered  in  this  work  and 
it  is  shown  how  the  model-based  approach  can  be  used  for  time  propagation  of  delayed 
measurements  to  obtain  estimates  of  current  positions  of  other  agents  in  the  network.  Performance 
bounds  on  the  tracking  error  have  been  obtained,  which  are  functions  of  the  communication 
topology,  the  event  thresholds,  and  the  time  delays. 


I.  INTRODUCTION. 

Cooperative  control  problems  have  received  significant  attention  in  the  last  years.  These 
problems  typically  require  a  group  of  mobile  vehicles,  also  called  agents,  to  perform  a 
coordinated  task  by  exchanging  relevant  infonnation  according  to  either  a  fixed  or  a  time- 
varying  communication  topology.  Most  of  the  results  concerning  agents  with  continuous  time 
dynamics  require  each  agent  to  continuously  transmit  measurements  to  its  neighbors  [1],  [2]. 
However,  in  practical  problems  there  exist  computation  and,  especially,  communication 
limitations  that  impose  constraints  into  how  frequently  an  agent  is  able  to  receive  measurement 
updates  from  its  neighbors. 

Recently,  typical  problems  in  cooperative  control  have  been  studied  from  a  sampled-data 
perspective  [3] -[9]  where  periodic  samples  of  continuous  time  output  signals  are  transmitted  over 
a  communication  channel  and  are  also  used  to  compute  sampled  control  inputs.  The  work  in  [3] 
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considers  a  sampled-data  consensus  approach  for  double-integrator  dynamics  and  for  undirected 
graphs.  The  authors  of  [4]  consider  sampled-data  consensus,  and  two  communication  cases: 
synchronous  and  asynchronous.  Reference  [5]  studied  the  consensus  problem  with  sampled 
information  and  time-varying  communication  delays.  Their  approach  considers  switching 
communication  topologies  and  consensus  is  obtained  for  small  enough  sampling  periods  and  if 
there  exists  a  frequent  enough  directed  spanning  tree  in  the  presence  of  delays. 

In  the  present  paper,  we  extend  the  results  on  sampled-data  cooperative  tracking  provided  in 
[6],  in  order  to  reduce  communication  rates  between  agents  and  to  improve  performance  as 
measured  by  the  tracking  errors.  We  make  use  of  the  Model-Based  Event-Triggered  (MB-ET) 
control  strategy  [10]-[1 1],  [  14]-[  1 5]  in  which  each  agent  implements  a  nominal  dynamical  model 
of  other  agents  in  the  network.  The  states  of  the  models  provide  an  estimate  of  the  real  positions 
of  neighbor  vehicles  and  these  estimated  variables  are  used  to  compute  the  control  inputs 
between  measurement  updates.  Each  agent  generates  its  corresponding  updates  based  on  local 
error  which  measure  the  difference  between  real  and  estimated  positions.  The  use  of  event- 
triggered  control  in  cooperative  control  has  been  used  by  several  authors  [12],  [13].  In  these 
papers  a  Zero-Order-Hold  (ZOH)  model  is  used,  that  is,  the  measurement  updates  received  by 
each  agent  are  held  constant  until  new  measurements  arrive.  In  the  approach  shown  in  this  paper, 
the  updates  received  from  neighbors  are  used  by  each  agent  to  update  an  internal  model  of  the 
neighbor  and  provide  an  estimate  of  the  neighbor’s  position  between  measurement  updates. 

The  paper  is  organized  as  follows:  Section  II  states  the  problem.  Section  III  presents  the  main 
results  of  this  note.  Examples  are  provided  in  Section  IV  and  conclusions  are  given  in  Section  V. 

II.  PROBLEM  STATEMENT. 

The  problem  studied  in  this  paper  is  the  dynamic  leader  coordinated  tracking  problem  with 
limited  communication  where  a  group  of  n  agents  receive  infonnation  about  the  leader’s  position 
at  some  time  instants  and  they  are  required  to  track  the  leader’s  position  as  close  as  possible.  The 
approach  we  follow  is  a  model-based  approach  in  which  each  agent  is  equipped  with  dynamical 
models  of  other  agents  in  the  group,  a  model  of  the  leader,  and  its  own  model  as  well. 

Consider  a  group  of  n  agents,  labeled  as  followers.  Assume  the  dynamics  of  each  agent  are 
given  by  the  following  continuous  time  dynamics: 
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ri(t)  =  ul(t),  i  =  1,2, 

The  control  inputs  are  given  in  sampled-data  form  as  follows: 

ui(t)  =  uj[k],  for  kT  <t  <(k  +  \)T.  (2) 

We  use  the  sampling  time  T>  0  to  discretize  the  continuous  time  dynamics  (1)  to  obtain: 

ri{k  +  \]  =  ri{k]  +  Tui[k].  (3) 

Note  that  the  sampling  time  T  is  only  used  to  obtain  a  discrete  time  equivalent  of  the  continuous 
time  agent  dynamics  but  it  is  not  used  to  establish  a  periodic  communication  between  agents.  At 
every  time  k  each  agent  will  estimate  the  positions  of  other  agents  and  its  own,  compute  its  local 
error,  and  update  its  local  input. 

The  discrete  time  models  that  will  be  implemented  by  each  agent  are  represented  by: 

r«)[k  +  l]  =  riU)[k\  +  Tu«)[k\,  /  =  0,1,2,...,/!  (4) 

where  f/  1  [A  J  represents  the  model  state  of  agent  i  estimated  by  agent  /  and  u(f[k\  is  the 
corresponding  model  control  input.  Note  that  a  model  of  agent  0  (the  leader)  is  also  implemented 
by  each  agent.  The  difference  between  the  inputs  ufk]  and  ufk]  will  become  apparent  in  the 


following  section. 

In  order  to  determine  the  time  instants  at  which  a  given  agent  needs  to  broadcast  its  current 
measured  position  we  use  an  event-triggered  strategy.  The  events  are  triggered  by  the  size  of  the 
local  agent  state  error.  The  local  state  errors  are  given  by: 


^[k\  =  ri[k\-rf[k\,  i  =  0,1,2 


(Or 


,n. 


(5) 


The  local  errors  measure  the  difference  between  the  real  agent  position  and  the  position 
estimated  by  its  own  local  model.  When  agent  i  decides  to  transmit  a  measurement  update  then 
we  have: 


?(>') 


(6) 


where  k 1  represents  the  update  instants  for  agent  i,  i.e.  the  time  instants  when  agent  i  updates  its 
model  and  transmits  this  update;  we  use  this  notation  to  emphasize  the  fact  that  the  updates  do 
not  take  place,  in  general,  at  every  time  k,  but  at  some  irregular  instants  k‘  . 

Each  agent  will  transmit  its  current  measured  position  to  the  rest  of  the  agents  if  its  local  error 
is  larger  than  a  predefined  threshold  otherwise  the  agent  will  not  attempt  to  send  measurements 
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since  the  current  estimate  is  close  to  its  real  position.  When  an  agent  updates  its  model  the  error 
is  equal  to  zero  at  that  time  instant  following  the  update  (6).  Then  the  following  holds: 

s\l)[k\<a  i  =  0,1,2, (7) 

III.  PROPORTIONAL-LIKE  CONTROLLER. 

In  the  MB-ET  framework  the  leader  implements  a  model  of  its  own  dynamics.  Suppose  that  the 
leader  dynamics  can  be  represented  by: 

r0[k  +  \\  =  r0[k]  +  Tu0[k]  (8) 

for  unknown  w0[k] .  Since  the  input  for  agent  0  is  completely  unknown,  a  reasonable  and  simple 

model  for  the  leader  is  a  ZOH  model,  that  is,  the  leader  will  hold  its  latest  transmitted 
measurement.  Note  that  the  leader  only  implements  its  own  model,  the  ZOH  model  based  on  its 
broadcasted  measurements.  Since  the  leader  does  not  receive  updates  from  any  other  agent  it 
does  not  need  to  implement  models  of  other  agents.  Under  these  circumstances,  the  local  state 
error  and  the  ZOH  model  corresponding  to  the  leader  are  given  by: 

=  *[*]-?»[*] 

=  forkl<k<kl, 

where  k°,ju  =  1,2,...  represents  the  update  instants  corresponding  to  agent  0,  that  is,  the  time 

instants  when  the  leader  transmits  a  measurement  and  updates  its  model. 

In  this  work  we  consider  a  communication  topology  that  can  be  represented  by  a  directed  tree 
graph  with  root  at  the  leader.  A  simple  example  of  a  group  of  agents  transmitting  infonnation 
using  this  type  of  graph  is  shown  in  Fig.  1 .  In  this  configuration  each  agent  i  in  the  network  does 
not  need  to  implement  models  of  every  other  agent,  but  only  of  those  in  the  directed  path  0-/.  If 
agent  i  needs  to  send  measurements  to  at  least  one  other  agent,  it  implements  models  of  all 
agents  on  this  directed  path  including  itself.  However,  if  agent  i  is  not  sending  measurements  to 
any  other  agent  (i.e.  agent  i  is  a  leaf),  the  agent  will  not  implement  a  model  of  itself.  For 
instance,  agent  3  in  Fig.  1  is  the  destination  node  of  the  directed  path  0-1-3,  then,  agent  3  needs 
to  implement  models  of  agents  0,  1,  and  3.  Agent  6  only  needs  to  implement  models  of  agents  0, 
1 ,  and  3  since  it  never  sends  updates  to  any  other  agent. 
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Fig.  1.  A  group  of  n  followers  tracking  the  position  of  a  leader  using  a  directed 
tree  communication  graph. 


The  control  input  for  the  z-th  agent  is  given  as  a  function  of  the  agent’s  own  position  r,  and  the 
estimates  of  the  positions  of  the  remaining  agents  as  given  by  its  own  models.  The  control  inputs 
Uj[k\  are  given  by: 

«i[^]  =  -2^(^-[^]-OcoM)-a/o(^OT-^o°OT)  (10) 

j= i 

where  <x  is  the  (i,j)  entry  of  the  adjacency  matrix  associated  with  the  follower’s  communication 

graph  and  aj0  >  0  if  the  leader  is  a  neighbor  of  follower  i  and  ai0  =  0  otherwise.  Each  agent 

computes  its  own  control  input  since  it  has  access  to  its  real  position  r,[k\. 

The  model  control  inputs  used  by  agent  /  for  each  one  of  the  agents  that  agent  /  is  receiving 
information  from  (also  including  its  own  model  control  input,  i.e.  l=i)  can  be  described  by: 

fm).  do 

j= i 

The  model  control  inputs  can  be  computed  by  all  agents  in  the  network  since  every  agent 
implements  the  models  (4)  of  other  agents.  Equations  (11)  are  only  a  function  of  model 
variables;  therefore,  they  can  be  computed  at  every  node. 
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A.  Time  propagation. 

The  model-based  approach  is  also  useful  for  time  propagation  of  delayed  measurements  [10].  For 
instance,  the  model  implemented  by  agent  3  in  Fig.  1  is  given  by: 

>o(3)[£  +  l]l  M  0  0  IpfW 

r0)[k  +  \]  =  f{3)[k  + 1]  =  Tal0  1-Tal0  0  ff\k]  =  A/y)[k]  (12) 

f3(3)[k  +  l]  L  0  Tan  1-FaJ  f3(3)[k] 

which  contains  all  models  that  agent  3  needs  to  implement,  i.e.  models  of  agents  0,  1,  and  3. 
Suppose  an  update  is  triggered  at  agent  1  at  time  k,  which  means  that  agent  1  updates  its  model 
variables  using  its  real  position  as  follows: 


r^[k]]Jr^[k] 

n0>[k]\  [rx[k]  ' 


(13) 


The  entire  model  state  fll[A]  is  transmitted  by  agent  1  at  time  k.  The  update  (13)  will  be 
received  by  agent  3  after  dn  sample  delays,  at  time  k  +  d3 , .  Instead  of  updating  the 
corresponding  part  of  its  model  using  the  delayed  measurements  f(ll[A:]  agent  3  estimates  the 
current  value  of  ril)[k  +  d31\  using  the  portion  of  its  model  (12)  corresponding  to  r0  and  rx ,  that 
is,  agent  3  computes  [r^][k  +  d3X\,  q(3)[k  +  , J J7  using  the  initial  conditions 
[r0(3)[£],  q(3)[A:]]r  =[f0(1)[A:],  rxu[kY[  and  using  the  following  state  equation: 


>0<3)[/c  +  l]l_r  1  0  irf0(3)[A:]" 

/i(3)[£  +  1]J  “  \_Tal0  l-Ta10\[r^[k]  ' 


(14) 


Since  (14)  involves  simple  matrix  operations  the  estimate  [r0(3)[£  +  dn\,  )\{'\k  +  c/3, ]J7  can  be 
obtained  instantaneously  or  with  negligible  delay  after  the  delayed  information  has  been  received 
by  agent  3.  Then,  the  following  update  takes  place: 

>o3)[£  +  «y]  r  r™\k  +  dj 

q(3)[k  +  diX\  =  kx°\k  +  dix\  (15) 

f™[k  +  d3l\ J  [r3[k  +  d3  J 

this  means  that  the  model  variables  of  agent  3  are  updated  using  the  estimates  of  current  model 
variables  of  agent  1  and  the  real  measurement  corresponding  to  its  current  local  position 
r3[k  +  <f31] .  This  update  also  generates  an  event  and  the  updated  model  variables  (15)  are 
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transmitted  to  agents  6  and  7  in  the  example  in  Fig.  1.  These  agents  perfonn  the  corresponding 
propagation  operations  when  they  receive  the  update  at  times  k  +  d3i  +  <7  63  and  k  +  dn+d73, 
respectively.  Agent  4  performs  similar  time  propagations  and  updates  as  agent  3  when  it  receives 
the  update  from  agent  1  at  time  k  +  d4l .  In  general,  the  delays  are  allowed  to  be  different  for  each 
communication  link. 

Note  that  transmission  events  are  generated  when  the  local  agent  error  reaches  its  threshold 
value  and  also  when  an  agent  receives  an  update  from  its  parent  such  as  agent  3  above;  agent  3 
received  an  update  from  agent  1  and  it  needs  to  resend  this  information  to  agents  6  and  7  since 
they  also  implement  a  model  of  agent  1 .  The  vector  of  infonnation  is  augmented  to  include  the 
current  position  of  agent  3  as  well.  Since  the  modeling  error  of  the  leader  is  the  only  uncertainty 
in  this  framework,  events  in  the  follower  nodes  are  generated  only  when  they  receive  updates 
from  parent  nodes. 

Assuming  that  the  delay  dtj  is  known  to  the  receiving  agent  i,  e.g.  by  time  stamping,  then  we 
have  that: 

r{i)[k  +  dij]  =  ?U)[k  +  dij ]  (16) 

that  is,  the  propagated  variables  at  node  i  are  equal  to  the  current  model  variables  at  node  j  after 
delay  di} .  This  is  expected  since  both,  model  j  and  the  propagation  equation  at  node  i,  use  the 

same  parameters;  then,  by  using  the  same  initial  conditions  at  time  k  both  of  them  will  produce 
the  same  response  if  the  delay  dtj  is  known  to  agent  i  (and  if  no  update  has  taken  place  at  node  j 

since  the  previous  update).  In  general,  and  especially  for  long  delays,  we  need  to  consider  the 
possible  updates  received  by  node  j,  for  the  worst  case  delay. 

B.  Main  Results. 

Define  the  tracking  error  for  each  follower  as  follows: 

Sm  =  rt[k\-r,[kl  (17) 

Let  ^[k]  =  [£j [k],...,^n[k]f  denote  the  tracking  errors  in  vector  fonn. 

Proposition  1.  The  tracking  error  vector  response  is  given  by: 

i=i  i= i 
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where  Q  =  (In -TL-T\Jdiag{aw,...,an0}),  A  =  [A0A\,  A0  =  [a10,...,an0]r ,  Aj  represents  the  matrix 
obtained  when  all  rows  in  A  are  replaced  with  zero  row  vectors  except  for  the  /th  row,  and 
X'p[k]  -  ( -r0[k  + 1]  +  r0[k]) ln .  A  and  L  are,  respectively,  the  adjacency  and  the  Laplacian  matrices 
associated  with  the  followers  communication  graph.  The  state  model  error  vectors  are  defined  by 
£°'>[k]  =  [£Q)[k],...,ejl,'>[k]]T .  Each  individual  error  £j)[k]  =  rj[k]-rj°[k],  represents  the 

difference  between  the  real  position  of  an  agent  j  compared  to  the  model  of  j  implemented  by 
agent  i. 

Proof.  The  tracking  error  dynamics  for  follower  i  can  be  described  by: 
f[k  +  Y\  =  r,[k  +  Y\-r0[k  +  \] 

=  r,[^]  +  r(-2^(r<[^]-r;,)m)-a<o(r<[^]-r0(0[^]))-r0[^  +  l]  (19) 

j= i 

=  -  (r0[A:  + 1]  -  r0[A:]). 

j= i 

The  tracking  errors  can  be  expressed  in  compact  fonn  as  follows: 
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(20) 


:(/„-rk-r^ag{a10,...,a,,J)^]-T(fi/1)[k]  +  fi2^[k]  +  ...  +  V",W)  +  ^;W- 


U2)r 


-('Or 


The  response  of  (20)  with  initial  conditions  2'[0J  and  with  inputs  £U)[k],  X’p[k]  can  be  directly 
described  as  in  (18). 

Theorem  2.  Assume  that  the  continuous  time  leader’s  position  r0(t)  satisfies  the  following: 


r0(t)\<r. 


(21) 
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For  any  7>0,  \r0(t  +  T)-r0(t)\  <  FT  holds,  which  is  equivalent  to  (when  considering  T  as  the 
discretization  sampling  interval) 


r0[k]-r0[k-l] 

T 


<  r 


(22) 


then,  the  maximum  tracking  error  of  the  n  followers  is  ultimately  bounded  by 


(s(i) +  r)\\(L  +  diag{al0,...,aii0}) 


-1 


(23) 


when  T  is  designed  to  satisfy 


T  <  mini/  . 

i=\,...,n  J 


(24) 


i= o 


Proof.  From  (18)  the  norm  of  the  tracking  error  can  be  expressed  as: 


IM.s  2‘Ml+r 


k- 1 


Zetv'r- i]+v!l[i-i]+-.+V'['-i]) 


1=0 


+  Tr 


k- 1 


/=0 


(25) 


Since  we  consider  a  directed  graph  every  follower  receives  updates  from  one  and  only  one  node. 
Then  we  only  need  to  consider  the  effect  of  the  estimation  errors  s+k]  for  which  at/  >  0  .  These 

errors  correspond  to  the  difference  between  neighbors  (agent  j  is  a  neighbor  of  agent  i). 

The  errors  are  given  by: 

sf[k]  =  r[k]  - +[k]  =  £{f[k]  +  +\k]  -  rf\k].  (26) 

For  the  first  tenn  in  the  right  hand  side  of  (26)  we  have  that  (7)  holds.  For  the  difference  between 
models  of  agent  i  and  its  neighbor  j  we  consider  the  response  of  the  entire  model  state.  If  an 
event  is  generated  by  agent  j  at  time  k  and  the  model  state  vector  rU)\k\  is  sent  to  agent  i,  we  are 


interested  in  the  largest  possible  difference  at  time  k  +  dtj  ,  where  dtj  represents  an  upper  bound 
on  the  communication  delay  from  agent  j  to  agent  i, 

(27) 


PJ '[k  +  df-P1 '[k  +  dy] 


r^fk  +  df-r^ik  +  df 


that  is,  the  difference  between  the  response  of  the  model  j  from  the  time  of  update  at  time  k  with 
initial  conditions  rij)\k\  and  the  time  propagation  of  the  measurement  r(j)[k]  by  agent  i,  where 
agent  i  assumes  that  no  other  event  has  been  generated  by  agent  j.  In  order  to  bound  the 
difference  (27)  we  consider  the  worst  case  scenario  where  the  leader  moves  using  its  maximum 
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velocity  and  (22)  holds  with  equality.  If  the  delay  is  small  such  that  the  rTdy  <  a  holds,  then  we 
have 


?0')r 


[k  +  diJ]-r^[k  +  dij] II  < 


Adiir  U) 
Aj  ra 


[k]-Aa/rU)[k] 


A,‘J[a  0...0]7 


(28) 


where  rai])[k]  =  r{J\k]  +  [cc  0 ...  0]r  and  A  represents  the  all  model  dynamics  used  by  agent  /  as 
illustrated  in  (12).  In  the  case  that  the  delays  are  larger  and  rTdy  >  a  holds,  then  additional 

events  have  been  generated  by  agent  j  before  the  measurement  r(r,[k]  reaches  agent  i.  For  this 
case  we  use  an  approximation  for  the  response  of  the  model  of  j  which  considers  the  greatest 
possible  difference  between  r^J)[k\  and  r0(o[£] ,  which  is  given  by  rTdy  +  a  ,  as  initial  condition 
of  agent  j  at  time  k.  An  illustration  of  these  responses  for  the  nodes  1  and  3  is  shown  in  Fig.  2. 


10  10.5  11  11.5  12  12.5  13  13.5  14  14.5 

Time  (seconds) 

Fig.  2.  Example  of  response  of  models  fU)[/:] ,  7'(1)[A']  of  adjacent  agents  1  and  3  with 
r=0.01  seconds,  constant  delay  of  3.3  seconds  and  worst  case  leader’s  velocity.  x<1)[A] 
represents  the  approximation  of  the  response  of  model  1 . 
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Figure  2  shows  the  approximation  of  rm[k\  =  [r^[k\  if\k]]T  denoted  by 
x(X)[k\  =  [xiQ)[k\  x^[k]]T  at  A=10  with  7=0.01  seconds,  a  constant  delay  of  3.3  seconds  (330 
sampling  times)  and  cc=  1. 

In  general,  considering  the  worst  case  leader’s  velocity,  we  have  that: 

|f0)[£  +  Zj-f(0[£  +  Zj  <«+  A^[7Td„+a  0...0f  (29) 

II  J  J  lloo  J  J  00 

for  any  two  agents  i  and  j,  where  j  is  a  neighbor  of  i.  Therefore,  a  bound  for  the  estimation  error 
(26)  at  any  time  k  is  given  by: 

ef[k]  <2 a+  A** [rTd  +a  0 ...  Of  (30) 

oo 

Define: 

s{,)  =  max  aib£^] [&]|  (31) 

i=\...n  1  '  I 

where  bt  denotes  the  neighbor  of  agent  i,  for  /=  1 . .  .n.  When  all  weights  are  equal  a  =  atj  then  (31) 
is  given  by: 

s{i)  =a  max  4° [it]  =  a(2a  +  A^rTd+a  0...0f  )  (32) 

i=i...n  1  \  J  oo  / 


where  d  =  max  h/J  is  the  largest  communication  delay  between  any  two  adjacent  agents. 

i=l...n  I  J  I 

By  considering  a  directed  graph  tree  we  can  write  (20)  in  the  following  form: 

&k+i\=Qm-Tm+x;[k\ 

where 


aw£^[k] 
a2b2Sb2  t^] 


(33) 


(34) 


(35) 
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n 


Since  the  leader  has  directed  paths  to  all  followers  and  0  <  T  <  min  1/  V  au  then,  by  lemma  8.3 


7=0 


in  [7],  Q  has  all  its  eigenvalues  within  the  unit  circle  and  lim  Qk  =  0 .  Additionally,  from  Lemma 


k  — ^co 


1.26  and  Lemma  1.28  in  [7],  we  have  that: 

llft-er'l  <36> 

and  the  maximum  tracking  error  of  the  n  followers,  ||c[/c  ]||  ,  is  ultimately  bounded  by  (23). 

The  state  error  for  the  leader  contains  significant  model  uncertainties  since  it  is  not  possible  to 
predict  the  behavior  for  this  particular  vehicle.  A  ZOH  model  is  implemented  by  the  followers 
and  also  by  the  leader  itself  and  it  was  given  in  (9).  The  state  error  for  the  leader  increases  in 
general  after  every  update  for  a  non-stationary  leader.  The  leader  is  able  to  store  its  latest 
broadcasted  update  and  compute  (9)  in  order  to  decide  when  a  new  update  should  be  broadcasted 
to  the  rest  of  the  agents.  Note  that  (7)  holds  for  any  agent  including  the  leader. 

Remark  1.  In  Theorem  2  and  similar  to  the  results  in  [6]-[7]  T  has  to  meet  certain  conditions 
for  stability,  i.e.  it  has  to  be  small  enough  so  the  overall  system  does  not  become  unstable.  Using 
the  MB-ET  framework  we  are  able  to  reduce  T  to  obtain  a  stable  matrix  Q  but  we  do  not  need  to 
communicate  every  T  seconds.  We  are  now  able  to  avoid  frequent  communication  by  choosing  a 
desired  error  threshold  a  while  maintaining  a  stable  overall  system. 

Remark  2.  An  additional  advantage  of  the  framework  presented  in  this  paper  compared  to  the 
sampled-data  approach  [6]-[7]  is  that  the  agents  are  not  required  to  transmit  measurements  to 
neighbors  at  the  same  time  instants,  i.e.  at  the  sampling  instants  kT,  k=  0,1,...  In  the  sampled-data 
approach  the  agents  are  not  required  to  communicate  continuously  but  all  of  them  use  the  same 
period  T  to  sample  and  transmit  their  measurements.  The  communication  scheme  in  the  present 
paper  allows  for  asynchronous  communication  in  the  sense  that  each  agent  determines  its  own 
transmitting  time  instants. 


IV.  EXAMPLES. 

Example  1.  Consider  a  group  of  n=l  vehicles  following  the  position  of  a  leader  (Agent  0)  and 
transmitting  information  according  to  the  graph  in  Fig.  1.  The  parameters  for  this  example  are  as 
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follows:  al0  =  1,  a10  =  8,  a31  =  3,  a4l  =  4,  a52  =  5,  a63 3  =  6,  a73  =  6.  The  sampling  time  is  IK). 01  sec 
and  the  threshold  is  a=l.  The  delay  bounds  for  each  communication  link  are  as  follows: 
d10  =  2.2sec,  d20  =  2.25 sec,  d3l  =  3sec,  d4X  =  3.05sec,  d52  =  3.1sec,  da  =  2.5sec,  d12  =  2.9sec. 

Fig.  3  shows  the  position  of  the  leader  and  followers  over  time  (top  row)  and  the  tracking 
error  (bottom  row)  for  a  leader  with  constant  velocity. 


0  10  20  30  40  50 

Time  (sec) 


0  10  20  30  40  50 

Time  (sec) 


Fig.  3.  Positions  and  tracking  errors  for  MB-ET  leader  follower  with  constant  leader  velocity. 


Example  2.  Consider  the  same  example  with  the  same  weights,  threshold  sampling  time,  and 
communication  delay  bounds.  Fig.  4  shows  the  position  of  the  vehicles  over  time  (top  row)  and 
the  tracking  errors  (bottom  row)  for  the  case  of  a  leader  with  time-varying  velocity. 


V.  CONCLUSIONS. 

The  results  in  this  paper  extend  previous  work  concerning  cooperative  control  and  tracking  of 
a  dynamic  leader.  Previous  results  required  a  small  enough  sampling  time  for  bounded  tracking 
error  otherwise  the  overall  system  becomes  unstable.  The  approach  in  this  paper  implements  an 
event-triggered  control  and  a  model-based  approach  that  provides  estimates  of  positions  of  other 
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agents.  The  overall  framework  is  able  to  reduce  transmission  of  information  among  agents  while 
maintaining  a  bounded  tracking  error.  Simulation  examples  show  the  effectiveness  of  this 
approach.  Future  work  will  consider  the  implementation  of  models  of  neighbors  only  (versus  the 
entire  directed  path)  to  reduce  model  order  and  complexity,  and  will  also  consider  more  general 
communication  topologies. 


Time  (sec) 


Time  (sec)  Time  (sec) 


Fig.  4.  Positions  and  tracking  errors  for  MB-ET  leader  follower  with  time-varying  leader  velocity. 
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